TECHNICAL 

LIBRARY 


lAQ-AOESniq  | 


TECHNICAL  REPORT  ARLCB-TR- 80009 


A  WIDE  RANGE  K  EXPRESSION  FOR  THE  C-SHAPED  SPECIMEN 


J.  KAPP 
J,  C,  Newman 
J .  H ,  Underwood 


March  1980 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BEN^T  WEAPONS  LABORATORY 
WATERVLIET ,  N.  Y.  12189 


AMCMS  No.  53970M63500 


PRON  No.  1A924154GGGG 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


183?6CIEDSj 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  author¬ 
ized  documents. 

The  use  of  trade  name(s)  and/or  manufacturer (s)  does  not  consti¬ 
tute  an  official  indorsement  or  approval. 


DISPOSITION 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 


to  the  originator. 


9  June  1980 


ERRATA  SHEET 
(Change  Notice) 

Cl  TO:  TECHNICAL  REPORT  ARLCB-TR-80009 


A  WIDE  RANGE  K  EXPRESSION  FOR  THE  C-SHAPED  SPECIMEN 


by 

J .  A .  Kapp 

J.  C.  Newman,  Jr. 

J.  H.  Underwood 


''Remove  the  Report  Documentation  Page,  DD 
from  the  above  subject  publication  dated 
and  insert  new  Report  Documentation  Page 


Form  1473, 
March  1980, 
inclosed . 


US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BENET  WEAPONS  LABORATORY 
WATERVLI ET ,  N.  Y.  12189 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

ARLCB-TR-80009 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

A  WIDE  RANGE  K  EXPRESSION  FOR  THE  C-SHAPED 

SPECIMEN 

S.  TYPE  OF  REPORT  &  PERIOD  COVERED 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORf  a) 

J.  A.  Kapp,  J.H.  Underwood 

J.C.  Newman  Jr.-  NASA  Langley  Research  Ctr. 
Hampton,  VA  23665 

8.  CONTRACT  OR  GRANT  NUMBER^ a) 

9  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Benet  Weapons  Laboratory 

Watervl'iet  Arsenal,  Watervliet,  NY  12189 
DRDAR-LCB-TL 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

AMCMS  No.  53970M63500 

PR0N  No.  IA924154GGGG 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  &  Development  Command 
Large  Caliber  Weapon  Systems  Laboratory 

Dover,  New  Jersey  07801 

12.  REPORT  DATE 

March  1980 

13.  NUMBER  OF  PAGES 

15 

14.  MONITORING  AGENCY  NAME  &  ADDRESS^//  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

ISa.  DECLASSI  F|  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

Presented  to  ASTM  Taskgroup  E24. 01.05,  30  Oct  79,  Pittsburgh,  PA. 

To  be  published  in  the  Journal  of  Testing  and  Evaluation. 

19.  KEY  WORDS  (Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

Calibration 

Fracture  Mechanics  Fracture  Testing 

Fracture  Properties  Toughness 

20.  AB5THACT  (Contfnus  an  reverse  sfxtim  ff  necessary  and  Identify  by  block  number) 

A  new  expression  has  been  developed  to  calculate  K  for  the  C-shaped  speci¬ 
men  over  a  wider  range  of  specimen  parameters,  namely  a/W,  X/W,  and  r^/r^, 
than  had  previously  been  available.  The  rationale  used  to  derive  the  x 
expression  was  to  utilize  known  stress  intensity  factor  solutions  for  short 
and  deep  cracks  to  develop  a  nondimens ional  form  of  K  which  approaches  finite 
values  as  a/W  goes  to  both  zero  and  one.  Numerical  K  solution  results  from 

prior  work  were  then  nondimensionalized  to  this  form  with  the  finite  limiting 

DD  t '22-73  M73  EDITION  OF  1  NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Deta  Entered) 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE (Wh an  Data  Entered) 


20.  ABSTRACT  (Cont'd) 

values,  and  the  dependence  of  K  on  a/W  was  determined  by  multi-variable, 
linear  regression.  The  final  expression  agrees  with  the  numerical 
K  solutions  within  ±  1.0%  for  .45  £  a/W  <  .55  for  all  r  /r  and  X/W  of 
either  0  or  .5;  withini:  1.5%  for  .2  <_  a/W  <_  1  for  all  ^  /r  z  and  x/W 
equal  to  0  or  .5;  and  within  about  +  3%,  for  .2  £_  a/W  £.  1  for  all  i^/r^ 
and  0  X/W  ^  1.  The  accuracy  of  this  expression  will  allow  expanded  use 
of  the  C-shaped  specimen  for  R-curve  determination  and  fatigue  crack  growth 
rate  testing. 


SECURITY  CLASSIFICATION  OF  THIS  PAGEfWTien  Data  Entared) 


CONTENTS 


Page 

INTRODUCTION  1 

PROCEDURE  1 

COMPARISON  WITH  COLLOCATION  DATA  5 

RESULTS  6 

REFERENCES  12 

TABLES 

1.  Comparison  of  f(a/w)  Values  from  Equation.  (10)  with  Corres-  9 

ponding  Values,  fc,  from  Collocation  and  Limit  Solutions. 

2.  Comparison  of  f(a/w)  Values  from  Equation.  (10)  with  Corres-  10 

ponding  Values,  fc,  from  Reference  4. 

3.  Values  of  f(a/w)  for  the  Range  of  a/w  used  in  K^c  Testing.  11 

4.  Values  of  f(a/w)  for  Wide  Range  of  a/w.  11 

ILLUSTRATION 

Fig.  1.  The  C-Shaped  Specimen,  Indicating  Geometric  Parameters  8 


i 


INTRODUCTION 

Recently  the  C-shaped  specimen  has  been  included  as  a  standard  specimen 
in  ASTM  Test  for  Plane-Strain  Fracture  Toughness  of  Metallic  Materials, 
E-399-78.  The  stress  intensity  factor  (K)  expression  for  the  C-shaped  spec¬ 
imen  in  the  ASTM  Test  was  obtained  by  Underwood  and  Kendall  [1]  by  fitting  a 
three  parameter  expression  to  boundary  collocation  results.  The  expression 
applies  over  the  following  ranges  of  the  three  parameters: 

0. 3  <_  a/W  <_  0. 7  ,  0_<X/W<_0.7  ,  0  <  <  l.°  , 

where  a  is  crack  length,  W  is  specimen  depth,  and  r^  and  r^  are  respectively 
the  inner  and  outer  radius  of  the  specimen  as  shown  in  Figure  1.  For  fatigue 
crack  growth  rate  testing  and  R-curve  measurements  using  the  C-shaped  specimen, 
K  must  be  known  over  a  wider  range  of  a/W  than  that  of  the  current  expression. 
The  purpose  of  this  paper  is  to  develop  a  wider  range  K  expression  for  the 
C-shaped  specimen. 

PROCEDURE 

Srawley  [2]  has  developed  wide  range  K  expressions  for  two  Kj^  specimens 
used  in  ASTM  Method  E-399-78,  the  compact  specimen  and  the  bend  specimen.  An 
approach  similar  to  that  used  by  Srawley  was  used  here  to  develop  a  K  expres¬ 
sion  for  the  C-shaped  specimen.  Since  K  varies  from  zero  to  infinity  as  a/W 
goes  from  zero  to  one  for  the  C-shaped  specimen,  it  is  very  difficult  to 
determine  an  expression  which  accurately  represents  the  variation  in  K  with 
crack  depth.  This  problem  can  be  circumvented  by  using  a  nondimensional 
form  of  K  which  has  the  correct  limiting  values  as  a/W  approaches  both  zero 
and  one.  Known  K  solutions  for  a  finite  crack  in  a  semi-infinite  plane  and 
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for  a  semi -infinite  crack  in  a  semi-infinite  plane  were  applied  as  limiting 
cases  for  short  and  deep  cracks  respectively.  This  resulted  in  an  expression 
which,  for  both  limits  of  a/W,  converges  to  the  correct  limiting  value.  The 
a/W  dependence  between  the  limits  was  described  by  fitting  a  polynomial  in 
a/W  to  the  available  K  solutions  [3-6]  for  the  C-shaped  specimen. 

Because  the  C-shaped  specimen  is  intended  for  use  with  a  range  of  hollow 
cylinder  geometries,  the  K  expression  must  account  for  load  eccentricity,  X/W, 
and  the  radius  ratio,  r^/r^,  in  addition  to  a/W.  The  X/W  dependence  was 
determined  by  examining  the  limiting  K  solutions,  and  the  r./r  dependence  was 
found  using  the  numerical  K  results. 

First  we  consider  the  deep  crack  limit  which  is  a  function  of  the  resul¬ 
tant  normal  force,  P,  and  the  resultant  moment,  M,  acting  on  the  uncracked 
ligament  of  the  specimen.  If  P  is  assumed  to  act  at  the  center  of  the  uncracked 
ligament  with  dimension  (W-a),  the  stress  intensity  factor  is  given  as  [7]: 

Kt  =  0.464  ■  -2P  .  (1] 

1  B /it  (W-a) 

For  the  case  of  the  resultant  bending  moment,  M,  K  is  [7]'. 


Kj  =  3.975 


M 


“3/2  • 


(2) 


B(W-a) ' 

The  total  stress  intensity  factor  as  the  crack  depth  approaches  through 
penetration  of  the  specimen  is  the  superposition  of  equations  (1)  and  (2) . 

By  calculating  M  in  terms  of  C-shaped  specimen  parameters  and  after  some 
algebra,  the  limiting  value  as  a/W  approaches  1,  of  the  familiar  nondimen- 
sional  form  of  K  is  determined. 
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(3) 


lim  KB/ftf  1.325 (3X/W  +  1.926  +  1.104a/W) 

a/W"1  P  Cl-a/W)3''2 

By  rearranging  equation  (3) ,  a  nondimensional  form  of  K  which  approaches  a 
finite  value  for  all  X/W  and  as  a/W  goes  to  1  is  obtained. 

3/2 


lim 


KBvW  (1-a/W) 


a/W  *  1  P (3X/W  +  1.926  +  1 . 104a/W) 


=  1.325 


(4) 


To  determine  a  nondimensional  form  of  stress  intensity  factor  that  has  a 
finite  limit  for  short  cracks  a  different  K  solution  must  be  used,  the  well 
known  Wigglesworth  solution, 

Kj  =  1.12  a/ira  (5) 

where  a  in  this  case  is  the  total  stress  acting  perpendicular  to  the  plane  of 
the  crack  at  r^  in  an  uncracked  C-shaped  specimen.  This  stress  has  two  com¬ 
ponents,  a  normal  component  and  a  bending  component,  of  which  the  bending 
stress  is  strongly  dependent  on  the  curvature  (r^/r^l  of  the  specimen.  Since 
the  deep  crack  limit  is  unaffected  by  r^/^,  the  effects  of  curvature  must- 
disappear  as  the  relative  crack  depth  increases  to  1.  The  r^/r.,  dependence 
cannot  be  totally  accounted  for  by  analysis  of  the  short  crack  limit.  There¬ 
fore,  the  stress  o,  used  in  equation  (5)  is  initially  calculated  assuming  no 
effects  of  curvature  and  the  r^/^  dependence  will  be  determined  separately. 

Using  linear  bending  theory  and  basic  mechanics  of  materials  concepts, 
the  stress  o  is  easily  determined  in  terms  of  C-shaped  specimen  parameters. 
Substituting  this  in  equation  (5),  the  limiting  value  of  the  familiar  non- 
dimensional  form  of  stress  intensity  factor  is 
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(6) 


lim 

a/W  -*•  0 


KBv¥  _  2.24/jf  >^a7W(3X/W  +  2  +  a/W) 
P  (1-a/W) 2 


Upon  rearranging  equation  (6)  a  nondimensional  form  of  K  is  obtained  that  has 
a  finite  limit  for  all  X/W  as  a/W  goes  to  zero. 


lim 


KB^d-a/W) 


=  2.24/it 


a/W  0  Pv^a7W(3X/W  +  2  +  a/W) 

Although  the  form  of  equations  (4)  and  (7)  are  somewhat  different,  some 
changes  can  be  made  to  both  equations  that  do  not  change  these  limits.  B&c  ause 
/a/W  approaches  1  as  a/W  approaches  1,  this  term  may  be  included  in  the  demon- 
inator  of  equation  (4)  without  disturbing  that  limit.  Also  for  a/W  close  to 
zero,  (1-a/W)2  =  (1-a/W)  ^  to  a  first  order  approximation,  thus  (1-a/W) 
can  replace  (1-a/W)2  in  the  numerator  of  equation  (7)  leaving  that  limit 
unchanged.  The  remaining  difference  between  the  two  equations  does  affect 
the  limits.  Since  the  form  of  equation  (4)  gives  a  better  description  of  the 
X/W  dependence  of  the  collocation  solutions  than  that  of  equation  (7),  it  was 
decided  to  use  a  simplified  form  of  equation  (4)  as  the  form  of  nondimensional 
K  used  to  fit  the  collocation  data.  Thus  the  final  form  of  the  wide  range 
expression  is: 

KBv¥(l-a/W)5/2 


(7) 


.1/2 


=  f (a/W)  *  g(ri/r2  ,  a/W) 


P(a/W)  (3X/W  +  1.9  +  l.la/W) 

The  functions  f  and  g  were  then  fit  to  the  collocation  data  with  the  condi- 


(8) 


tions  that  for  the  limits  a/W  1  and  a/W  0  the  product  of  f  and  g  equals 
the  values  in  equations  (4)  and  (7)  respectively. 


4 


The  function  g  was  difficult  to  determine  exactly.  For  all  r  /^  as 
a/W  approaches  1,  g  must  approach  1,  and  as  a/W  approaches  0,  g  is  the  ratio 
of  the  stress  normal  to  the  plane  of  the  crack  at  r^  using  curved  beam  theory 
to  that  stress  acting  at  using  straight  beam  theory.  This  ratio  was  deter¬ 
mined  and  found  to  be  highly  complex  and  has  a  strong  dependence  on  X/W  in 
addition  to  r^/r^.  Therefore,  g  was  developed  by  fitting  a  polynomial  that 
describes  the  r^/r^  ^ePenc^ence  the  available  collocation  results  for  the 
standard  C-shaped  specimens,  that  is  for  X/W  =  0  and  0.5.  For  a/W  ranging 
from  .2  to  1,  g  was  found  to  be: 

gOyr^a/W)  =  1  +  . 25 (1-a/W) 2  (1-r^/r^)  .  (9) 

The  function  f  was  determined  from  the  available  collocation  data  and  the 
deep  crack  limit  for  the  standard  C-shaped  specimens.  Using  multivariable 
linear  regression,  these  data  were  fit  to  the  following  cubic  polynomial,  for 
a/W  from  0.2  to  1. 

f (a/W)  =  3.74  -  6.30  a/W  +  6.32(a/W)2  -  2.43(a/W)3  (10) 

COMPARISON  WITH  COLLOCATION  DATA 

To  illustrate  the  goodness  of  fit  of  equation  (10),  values  of  f(a/W) 

calculated  from  the  collocation  results  for  standard  specimens  are  compared 

in  Table  1  to  the  corresponding  values  determined  from  equation  (10) .  For 

comparison  we  have  defined  an  error  function  E,  which  is  the  difference 

between  the  value  of  f(a/W)  from  equation  (10)  and  the  value  f^  from  the 

collocation  results  divided  by  the  value  of  f  . 

c 
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Based  on  the  comparisons  in  Table  1  and  on  the  information  from  the 
limit  solutions,  equation  (10)  agrees  with  all  these  solutions  within  ±  1.5% 
for  .2  <  a/W  <_  1  and  within  ±  1%  for  .45  <_  a/W  <  .55  for  all  r  /r  and  X/W  of 
either  0  or  .5. 

Comparisons  are  made  in  Table  2  between  equation  (10)  and  other  collo¬ 
cation  results  for  nonstandard  specimens.  Analysis  of  the  relative  errors 

in  these  cases  indicate  that  the  derived  K  expression  should  not  be  used  for 
values 

specimens  with  X/W^greater  than  about  1.  As  X/W  beamnes  largo,  the  banding 
component  of  the  stress  distribution  becomes  dominant.  Since  the  bending 
component  is  significantly  affected  by  the  curvature  of  the  specimen,  the 
variation  in  K  could  have  been  explained  by  the  function  g  derived  above. 
However,  g  was  not  obtained  exactly,  but  by  fitting  to  the  numerical  solu¬ 
tions  for  the  standard  specimens,  so  good  agreement  can  not  be  expected  for 
X/W  larger  than  0.5. 

RESULTS 

To  calculate  K  using  equations  (8),  (9),  and  (10)  is  cumbersome.  To 
simplfy  K  determination,  equations  (8),  (9),  and  (10)  can  be  arranged  to  a 
more  convenient  form: 

K  =  -j-L  [3  X/W  +  1.9  +  1.1  a/W][l  +  .  25  (1-a/W)  2  (l-r^)  ]  F  (a/W)  ,  (11) 

where 

/a/w 

F  (a/W)  =  ---"3/2  C3*74  -  6.30  a/W  +  6.32  (a/W) 2  -  2. 43  (a/W)3) 

Numerical  values  of  F(a/W)  have  been  determined  for  two  ranges  of  a/W  and 
are  presented  in  Tables  3  and  4. 
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Equation  (11)  agrees  with  the  numerical  solutions  available  [3-6]  within 
±  1.0%for  .45  a/W  .55  for  all  r^/r 2  and  X/W  of  either  0  or  .5,  within 

±  1.5%  for  .2  a/W  <_  1  for  all  r  /r  anc*  x/w  equal  to  0  or  .5;  and  within 

about  ±  3%  for  .2  <_  a/W  1  for  all  r^/r2  anc*  ®  i  ^ •  The  x  expression 

of  equation  (11)  is  as  accurate  as  that  which  is  currently  given  in  ASTM 
Method  E-399-78.  The  new  expression  can  be  used  for  other  fracture  testing 
over  a  wider  range  of  a/W  than  that  of  the  current  expression. 
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Fig.  1,  The  C-Shaped  Specimen,  Indicating  Geometric  Parameters 


8 


TABLE  1.  COMPARISON  OF  f(a/W)  VALUES  FROM  EQ.  (10)  WITH  CORRESPONDING  VALUES,  f„, 

FROM  COLLOCATION  AND  LIMIT  SOLUTIONS 
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TABLE  2.  COMPARISON  OF  f  (a/W)  VALUES  FROM  EQ.  (10)  WITH 
CORRESPONDING  VALUES,  f  ,  FROM  REFERENCE  4 
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TABLE  3. 

VALUES  OF  F (a/W) 

FOR  THE  RANGE  OF  a/W  USED 

IN  Klc  TESTING 

a/W 

F (a/W) 

a/W 

F (a/W) 

.450 

3.23 

.500 

3.73 

.455 

3.27 

.505 

3.79 

.460 

3.32 

.510 

3.85 

.465 

3.37 

.515 

3.91 

.470 

3.42 

.520 

3.97 

.475 

3.47 

.525 

4.03 

.480 

3.52 

.530 

4.10 

.485 

3.57 

.535 

4.17 

.490 

'3.62 

.540 

4.24 

.495 

3.68 

.545 

4.31 

.550 

4.38 

TABLE  4.  VALUES  OF  F(a/W) 

FOR  WIDE  RANGE  OF  a/W 

a/W 

F (a/W) 

a/W 

F (a/W) 

.20 

1.70 

.60 

5.24 

.25 

1.94 

.65 

6.42 

.30 

2.20 

.70 

8.11 

.35 

2.49 

.75 

10.70 

.40 

2.82 

.80 

15.01 

.45 

3.23 

.85 

23.15 

.50 

3.73 

.90 

42.53 

.55 

4.38 

.95 

119.90 
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